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Orbital Kondo effect in a system of two single-level quantum dots attached to external electron 
reservoirs is considered theoretically. The dots are coupled via direct hoping term and Coulomb 
interaction. The Kondo temperature is evaluated from the scaling approach and slave boson tech- 
nique. The later method is also used to calculate linear conductance of the system. Nonlinear 
conductance, in turn, is calculated in terms of the nonequilibrium Green function formalism. 
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I. INTRODUCTION 

Kondo effect in electronic transport through quantum 
dots (QDs) strongly coupled to external leads is a many 
body phenomenon which has been extensively studied in 
the last two decades 1 - - — . Spin fluctuations in the dot, 
generated by coupling of the dot to external leads, give 
rise to a narrow peak in the dot's density of states (DOS) 
at the Fermi level. This Rado-Suhl resonance results in 
enhanced transmission through the dot, and leads to the 
unitary limit of the linear conductance G at zero tem- 
perature, G = 2e 2 /h. The enhanced transmission is sup- 
pressed when a bias voltage is applied, and this leads to 
the so-called zero-bias (Kondo) anomaly in differential 
conductance. The above described phenomenon arises 
from the two-fold spin degeneracy, and is often referred 
to as the spin Kondo effect. However, the Kondo phe- 
nomenon may also appear when the spin degree of free- 
dom is replaced by any two- valued quantum number, e.g. 
the one associated with an orbital degree of freedom (the 
orbital Kondo effect) 1 -. A minimal realization of the or- 
bital (spinless) Kondo phenomenon requires two orbital 
discrete levels coupled to external lead s 13 ' 14 . This can 
be realized for instance in two single-level quantum dots 
coupled to external electrodes^— . Coherent superpo- 
sition of virtual tunneling events, in which one electron 
tunnels from the dot QD1 (QD2) to one of the leads and 
then simultaneously another electron tunnels to the dot 
QD2 (QD1), leads to the Kondo resonance at low tem- 
peratures. 

In this paper we consider theoretically the Kondo phe- 
nomenon in electronic transport through two QDs cou- 
pled, in general, via both Coulomb interaction and hop- 
ping term. To evaluate the level renormalization and 
Kondo temperature of the system we use the scaling ap- 
proach. The Kondo temperature is also evaluated from 
the slave boson technique. Additionally, the latter tech- 
nique is used to calculate the linear conductance. Then, 
the nonequilibrium Green function formalism is used to 
calculate the local density of states (LDOS) for both dots 
and transport characteristics (differential conductance) 
in the nonlinear response regimes. To calculate the rele- 
vant Green's functions from the corresponding equations 
of motion we apply the decoupling scheme introduced in 



Ref. [|. 

The orbital Kondo effect in double quantum dot 
(DQD) systems was analyzed e.g. in Ref. [l8|]. How- 
ever, our results arc different, and the key difference con- 
sists in a different symmetry of the couplings to exter- 
nal leads. Moreover, we use various techniques including 
scaling, slave boson, and nonequilibrium Green function 
formalisms. Apart from this, we apply a different method 
to evaluate the lesser Green function. The paper is or- 
ganized as follows. In section 2 we describe the model of 
a double quantum dot system. Renormalization of the 
dots' levels and the Kondo temperature are discussed in 
section 3 in terms of the scaling approach. The slave bo- 
son technique is briefly described in section 4 and is used 
there to estimate the Kondo temperature and calculate 
the linear conductance. Basic formula for nonequilibrium 
Green functions and the corresponding numerical results 
on the conductance and LDOS are presented and dis- 
cussed in section 5. Summary and final conclusions are 
given in section 6. 

II. MODEL 

We consider two coupled single-level quantum dots 
connected to nonmagnetic electron reservoirs as shown 
schematically in FigJTJ Each dot is attached to sepa- 
rate source and drain leads, so the tunneling paths via 
the two orbitals can be analyzed separately as in recent 
experiment er 15 ' 16 . We consider the case when each dot 
is coupled symmetrically to the leads, while the corre- 
sponding coupling strengths for both dots may be dif- 
ferent. Moreover, our considerations are limited to the 
case of spinless electrons, which can be realized exper- 
imentally for instance by applying a sufficiently strong 
external magnetic field lifting the spin degeneracy. 

The system under consideration can be described by 
the extended Anderson Hamiltonian of the general form 

H = Hlcads + HuQD + ^tunnel- (1) 

The first term, H\ ca ,d s , describes here the four leads in the 
non-interacting quasi-particle approximation, ifioads = 
Hli + Ht,2 + Hri + Hr2, with Hpi being the Hamiltonian 
of the left (/3 — L) and right (/3 = R) lead attached to the 
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FIG. 1: (color online) Schematic picture of the double quan- 
tum dot system coupled to external leads. 



ith dot (i = 1,2), Hffi = Y,k e ^i c kpi c ^ffi ( for ^ = L,R 
and i — 1,2). Here, c\ t . j3i (ckp%) is the creation (annihi- 
lation) operator of an electron with the wave vector k 
in the lead j3i, whereas ekpi denotes the corresponding 
single-particle energy. 

The second term of the Hamiltonian (1) describes the 
double quantum dot system, 



H 



DQD 



eid]di+t{d[d 2 +h.c.) + Unin2, (2) 



where n, = d\d. t is the particle number operator (i = 
1,2), €i is the discrete energy level of the i-th dot, t de- 
notes the inter-dot hopping parameter (assumed real), 
and U is the inter-dot Coulomb integral. 

The last term, Ht, of Hamiltonian (1) describes elec- 
tron tunneling between the leads and dots, and takes the 
form 



V d 



h.c. 



(3) 



where are the relevant tunneling matrix elements. 
Coupling of the dots to external leads can be param- 
eterized in terms of rf(e) = 2tt J2k V ik V ik*^( e ~ £ k^)- 
We assume that is constant within the electron band, 



r H e ) = r f = const for 6 e (~D,D), and (e) = oth- 
erwise. Here, 2D denotes the electron band width. We 
assume the dots are symmetrically coupled to the leads, 
rf = rf = aT, and rf = = T. The parameter a 
takes into account difference in the coupling of the two 
dots to external leads. Note, that for these parameters, 
each dot separately is coupled symmetrically to the two 
leads. 



III. LEVEL RENORMALIZATION AND KONDO 
TEMPERATURE 

Coupling of the dots to external leads gives rise to 
renormalization of the energy levels of both dots. In 
this section we use the scaling approach to derive some 
general formula for the renormalized levels in unbiased 



system. From the scaling equations we also estimate the 
relevant Kondo temperature. The derived results will 
be used subsequently for interpretation of the numerical 
results on electronic transport and LDOS. 



A. Renormalization of the QDs' levels 

Now, we apply the scaling technique to derive renor- 
malized dots' energy levels and begin with the limit of 
t = 0. In the scaling approach, the high-energy excited 
states (in the energy region of width 3D at the band 
edges) are removed, but their impact on the system is 
taken into account via renormalized parameters of the 
Hamiltonian. Here, we consider only second order pro- 
cesses, where the leads' electrons are scattered to the 
band edges and back. To perform scaling we assume 
f-i + U 3> D 3> \ After integrating out the band edge 
states we arrive at the following renormalized parame- 
ters; 
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TjSD 
2tt D ' 



(4) 



where Eq is the energy of empty DQD system (initially 
E = 0) and the index i = 1 for i — 2 and i = 2 for i = 1. 
Here, 1^ is defined as T t = Tf'+Tf. This procedure leads 
to the following scaling equation: 



dc, 



d\nD 

and to the level separation, 



2tt' 



Ae 
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ei - £2 = ei - (2 H r m — 

2vr \D 



(5) 



(6) 



where D is the band width at the end of scaling proce- 
dure. When ei = 62 = eo, the above equation shows that 
the initial degeneracy is generally lifted. 

Scaling in the presence of direct tunneling between the 
dots, t 7^ 0, is more complex in a general case. However, 
we restrict our considerations to some limiting cases, i.e. 
when the hoping term is weak, \t\/T <C 1, and when 
\t\/T ^> 1. If the bare dots' levels are degenerate, the di- 
rect hopping term generally lifts the degeneracy. The two 
eigenstates of the coupled quantum dots isolated from the 
leads correspond to the antibonding and bonding states, 
with the eigenenergies e± = (ei + £2)/2 ± -\/Ae 2 + 1 2 , 
where Ae = (ei — (.2)12. When the hoping term is small, 
one can first perform scaling of the bare dots' levels and 
then incorporate nonzero t by substituting ei(2) by £112) 
in the above expression for e±. The situation changes 
when tunneling coupling between the dots is larger than 
the dot-lead coupling. To find the relevant energy lev- 
els involved in the Kondo effect, one has to diagonalizc 
the dot's Hamiltonian first (transformation to the bond- 
ing and antibonding states), and then perform scaling for 
the energy levels e+ and e_. The corresponding scaling 
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equation has the form 



1.0 



de± 
d\nD 



2tt' 



(7) 



which is similar to Eq.©. However, the effective cou- 
pling of the new states to the leads acquires now the form 
r± = Ti + T2- Thus, the level separation in the limit of 
strong hoping term is independent of the couplings T± 
as they are both the same, T + = T_. 



B. Kondo temperature 

Now we evaluate the Kondo temperature for t — and 
(.\ = (.2 = £o using the 'poor man' scaling approach 26 . 
To pursue this method one has to derive first the Kondo 
Hamiltonian by performing the Schrieffer-Wolf transfor- 
mation. Then, the band width is reduced by eliminating 
states with energy D — 5D < |ek/3i| < D and introducing 
a new effective Kondo Hamiltonian, which has the same 
form as the initial one, but with renormalized parameters 
J+, J-, J z i and J Z 2- All information on the high energy 
excitations is incorporated into these renormalized pa- 
rameters. 

To apply the renormalization group procedure we first 
reformulate the definitions of the coupling parameters 
in the following way: Ti = 2oT = T(l — p) and T2 — 
2T= f(l+p), so that f = (a+l)r andp = (l-a)/(l + 
a) . After performing scaling procedure one arrives at the 
following scaling equations; 



d(pJ±) 
d\nD 
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= -pJ±(pJzl + pJz2), 



(8) 
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FIG. 2: Normalized Kondo temperature as a function of the 
parameter a, obtained from the scaling method for U = 5GT 
and from the slave boson technique for U — 00. The other 
parameters are eo = — 3.5F and t = 0. 



where instead of p we have spin polarization p of the 
leads. 

Variation of the Kondo temperature with the parame- 
ter a is shown in Fig|2] (dashed line). Tk reaches max- 
imum value for a = 1 (p = 0) and vanishes for a — > 
(p = 1). This behavior is similar to that for spin Kondo 
phenomenon in a QD coupled to ferromagnetic leads. We 
also note that our results are in agreement with those 
obtained in Ref.[28j], where the authors mapped spinless 
DQD system onto a spinful generalized Anderson model. 



d(pJzi) 
d\nD 



= -2(pJ±) 2 



(9) 



where p = pi = p 2 , with pi =J2p Ppu an d P/3i being the 
density of states in the lead j3i (f3 — L,R and i — 1, 2). To 
solve these equations we first find the scaling trajectories 
{pj±) 2 ~ {pJz\){pJzi) = and pJ z i - pJ z2 = const = 
pJzi ~ pJz2 — PP{Jzi + Jz2)- This allows us to write only 
one scaling equation instead of the two coupled equations 
© and ©, 



d(pJzi) 
d\nD 



= -2(pJ zi )[pJ zi T PP(J°zi + J°2)l 



(10) 



for i = 1,2. One actually continues the scaling process 
until D « IzbTk- Solving Eq. (|T0| one finds the Kondo 
temperature as the relevant scaling invariant, 



T K = Dexpl- 



1 arctanh(p) 1 



with p(J° zl + J« 2 ) = f M & +vy The above formula 
resembles the corresponding one for the Kondo temper- 
ature in a single QD coupled to ferromagnetic leads 2 ^, 



IV. SLAVE BOSON APPROACH 

To estimate the Kondo temperature and calculate con- 
ductance in the linear response regime, we apply now the 
slave boson technique for U — > oc 29 . This method relies 
on introducing auxiliary operators for the dots, and re- 
placing the electron creation and annihilation operators 
by fl b and b^fi, respectively. Here, b^ creates an empty 
state, whereas fj creates a singly occupied state with 
an electron in the «-th dot. To eliminate non-physical 
states, the following constraint has to be imposed on the 
new quasi-particles, 



(12) 



The above constraint prevents double occupancy of the 
system (the DQD system is either empty or singly occu- 
pied) . 

In the mean field approximation (MFA), the boson 
field b is replaced by an independent of time real num- 
ber, b{t) — > (b(t)) = b. This approximation, however, re- 
stricts considerations to the low bias regime (eV <C |e,-|). 
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Introducing now the following renormalized parameters: 

t = tb 2 , V"f k = Vfeb, and ej = e; + A, where A is the corre- 
sponding Lagrange multiplier, one can write the effective 
MF Hamiltonian as 

H MF = £ £ wU<** + £ e^/i + (f/t/ 2 + h.c.) 

k /5j t 

+ EE^ c U + h - c -) + a (& 2 - 1) . (i3) 

k /3i 

The unknown parameters, b and A, have to be found self- 
consistently from the following equations; 



(14) 



2^ 



(e-#A+^=0, (15) 
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FIG. 3: Linear conductance vs. dots' level position, ei = 
£2 = eo, for indicated values of a and t = 0, obtained from 
the slave boson method for U —t oo. The level splitting due 
to renormalization is tuned out by external gates. 



where ((fi\fj))f is the Fourier transform of the lesser 
Green function defined as G<{t,t') = ({fi(t)\fj (t')}}< = 
i(fi(t')fi(t)}. These equations follow from the constraint 
imposed on the slave boson field, Eq. (fl~2|) . and from the 
equation of motion for the slave boson operator. The 
lesser Green functions {{fi\f}))f as well as the retarded 
Green functions ((fi\fj)) r (the latter ones are required in 
the further calculations, too) have been determined from 
the corresponding equations of motion. 

The Kondo temperature can be introduced asii 



T, 



K 



(16) 



with r = & 2 (ri + T2) and ei = I2 — eo- Variation of the 
Kondo temperature (evaluated from the above equation) 
with the parameter a is shown in FigJ5] (solid line). 

To study charge transport we assume the same elec- 
trochemical potentials for the left leads and also equal 
electrochemical potentials of the right leads. The linear 
conductance is then calculated from the Landauer for- 
mula, in which the transmission matrix is taken at the 
Fermi level. More specifically linear conductance is given 
by the formula 



G 



,. dJ 



(17) 



where J is current calculated at temperature T = OK 32 . 
The slave boson technique in the form presented above, 
however, does not take into account the level renormal- 
ization described in the preceding sections. Therefore, to 
include this renormalization we replace the bare dot lev- 
els by the renormalized ones (keeping the notation used 
for the bare dot levels). Alternatively one may say that 
the renormalization is tuned out by external gate volt- 
ages. In Fig |3] the linear conductance is shown as a func- 
tion of the dots' energy level, e\ — £2 = £0i and for 



indicated values of the parameter a. The linear conduc- 
tance reaches the unitary limit for eo -C — F. This limit is 
achieved owing to the tuning out the level splitting due 
to renormalization. From this figure also follows, that 
the Kondo temperature decreases with decreasing a, in 
agreement with the above discussion and Fig. 2. 



NON-EQUILIBRIUM GREEN FUNCTION 
APPROACH 



Electric current flowing through a biased system is 
determined by nonequillibrium retarded, advanced, and 
lesser Green functions of the dots, and can be calculated 
from the formula derived by Meir et al2£. In turn, to cal- 
culate the retarded (advanced) Green functions G[j°^ (e) , 
we have applied the equation of motion method (EOM). 
Within this method one writes first the equation of mo- 
tion for the causal Green function (e) , which generates 
new Green functions. Then, one writes the equations of 
motion for these new Green functions, which in turn con- 
tain new higher-order Green functions. The latter ones 
have to be calculated approximately. To close the set of 
equations for the Green functions we have applied the 
decoupling scheme introduced in Ref.[4]. Although such 
an approximation does not describe properly the zero 
temperature limit, it is sufficient to describe the Kondo 
phenomenon close to the Kondo temperature. Detailed 
expressions for these Green functions are shown in the 
Appendix. 

The retarded/advanced Green functions contain occu- 
pation numbers, Uj, and the interdot correlators, = 
(d\di}, which can be calculated from the identities 
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(18) 



2ir 11 



(19) 



Thus, we still need the lesser Green functions G^(e). 
However, one can note that instead of Gfj(e), only 
/ deGfj(e) is needed. This quantity can be found ex- 
actly (in contrast to the approach based on the Ng's 
approximation^) and to do this we apply the Heisenberg 
equation of motion for the operators dj(t)di(t). Then, 
one takes average from the obtained equation and makes 
use of the fact that ((d/dt)dj(t)di(t)) = in a steady 
state. As a result one obtains the following equations; 



iYim 



= E/ ^f^i{Gl-Gl), 



(20) 



-{T l n fl +T inu ) 



t(rii - nj) + (q - e»)n 

= E/^( £ )[ r ^- r f G ^ 



(21) 



These equations, together with the appropriate equations 
for the retarded/advanced Green functions, have to be 
solved numerically in a self-consistent way. 

The basic transport characteristics of the system, like 
conductance and differential conductance can be calcu- 
lated numerically using the formulas derived above. The 
local density of states (LDOS) for the i-th dot can be 
calculated as 
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FIG. 4: Local density of states for the dots QD1 and QD2, 
calculated as a function of energy for indicated values of the 
parameter a. The other parameters are eo = — 3.5F, U = 50F 
and t = 0. 



D i = --Q[G r u {e)] 



(22) 



Numerical results for t = 



where $s[A] denotes the imaginary part of A. 

The approximation scheme used to calculate the 
nonequillibrium Green functions does not take into ac- 
count the level renormalization described in the preced- 
ing sections. Therefore, to take this renormalization into 
account we replace the bare dot levels by the renormal- 
ized ones (keeping the notation used for the bare dot 
levels), similarly as in the case of slave boson technique. 
However, one should bear in mind that the presented 
EOM approach just renormalizes the bare dot's energy 
levels due to real part of the corresponding self-energies. 
This renormalization can be seen looking at the position 
of the broad maximum in LDOS (see Fig|4]). However, 
the used decoupling scheme does not take it properly and 
thus does not lead to the expected splitting of the zero 
bias anomaly. In the following numerical calculations we 
assume equal dot energy levels, = e (for i = 1, 2) (e 
is measured from the Fermi level of the leads in equi- 
librium, fj,n = [iRi = 0). Apart from this, we assume 
e = -3.5T, the bandwidth 2D = 500r, and U = 50T. 



Let us start with the case when the dots are capac- 
itively coupled only, t = 0. The LDOS for both dots 
is plotted in FigJH The spectrum of each dot reveals 
two resonances corresponding to the dot level and its 
Coulomb counterpart (the latter not shown). Apart from 
this, a narrow peak emerges in the spectrum of each dot 
at the Fermi level of the leads. The intensity and width 
of this peak strongly depends on temperature, revealing 
all characteristic features typical of the Kondo resonance. 

The resonance in LDOS originates from the many body 
processes which occur in the low temperature regime. 
Since the conditions eo < fipt and fipi < eo+U are obeyed 
for the parameters assumed (Coulomb blockade regime), 
only a single electron can occupy the DQD system and 
sequential tunnelling processes are blocked. However, 
higher-order tunnelling events are still allowed. Let us 
assume that an electron initially occupies the dot QD1, 
and the system is in the Coulomb blockade regime. Due 
to the uncertainty principle, the electron from the dot 
QD1 can tunnel onto the Fermi level of one of the leads 
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attached to QD1, while an electron from the Fermi level 
of one of the leads attached to QD2 can tunnel to the 
dot QD2 in the time ?i/|eo|. Interference of many such 
events gives rise to the narrow peaks in LDOS at the 
Fermi level. 

For the fully symmetric model (a = 1), LDOS for the 
dot QD1 is the same as that for QD2. The situation 
is different for a ^ 1. As a decreases, the intensity of 
the Kondo peak in LDOS of the dot QD2 also decreases 
and disappears when a tends to zero. The opposite sit- 
uation occurs in the LDOS of the dot QD1, where the 
Kondo peak becomes more and more pronounced with 
decreasing a. This behavior is due to the fact that the 
intensity of the Kondo peak in LDOS of the dot QD1 
is mainly determined by the coupling strength between 
QD2 and the leads, while the Kondo peak for the dot 
QD2 is predominantly determined by coupling of the dot 
QD1 to the leads. Accordingly, the Kondo peak in LDOS 
of the dot QD1 (QD2) increases (decreases) with decreas- 
ing a, while the Kondo peaks of both dots are equal 
for a — 1. This can be understood taking into account 
the maximal value of the conductance per each channel 
(each dot), which is equal to G™- ax — e 2 /h. Estimating 
Qmax ^ Yn(E F )pi(E F ) at the Fermi level, one obtains 
pi{Ep) / pi{Ep) — 1/a, which explains the above behav- 
ior of the LDOS for QD1 and QD2. This behavior of the 
Kondo peaks in LDOS of both dots is similar to that in 
the case of spin Kondo effect, where each spin channel 
is coupled differently to the leads (when the leads are 
ferromagnetic) . 

For a <C 1, the Kondo peak for the dot QD1 becomes 
strongly asymmetric and the LDOS is totally suppressed 
for energies above the Fermi level, where the spectral 
function is equal to zero. This situation is similar to 
that reported in Ref.[l9|. Apart from this, position of 
the Kondo peak for QD2 slightly moves away from the 
Fermi level with increasing a (towards positive energies), 
and becomes asymmetric for all values of a. 

It is also worth to note that position of the broad max- 
imum (associated with the dot's level) in LDOS of QD1 
(the dot whose coupling to the leads changes with a) is al- 
most unchanged with tunning a, whereas position of the 
broad maximum in the LDOS of the dot QD2 (coupled to 
the leads with constant strength) varies with the parame- 
ter a. This becomes clear when considering the formulas 
for renormalized dots' energy levels, Eq.©. Apart from 
this, the intensity of the broad peak for the dot QD1 
decreases monotonically with increasing a, whereas the 
intensity of the broad peak in the LDOS of the dot QD2 
depends on a in a more complex way. When a — > 0, the 
broad peak in the LDOS of the dot QD1 is then most pro- 
nounced, whereas its Coulomb counterpart (not shown) 
is totally localized at e + U. 

As mentioned before, the resonances in LDOS lead to 
zero bias anomaly in the differential conductance of the 
DQD system. Here, nonlinear conductance is defined in 
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FIG. 5: Differential conductance for the dots QD1 and QD2, 
and the total differential conductance calculated for indicated 
values of a. The other parameters as in Fig[4] 



the following way: 



where J is current given by Meir and Wingreen 
formula^. This quantity is very important from prac- 
tical point of view as it is usually measured in QDs' ex- 
periments (to obtain basic transport properties of these 
systems) 1 . In Fig.© we show the differential conduc- 
tance of both dots as a function of the bias voltage. For 
a fully symmetric system, the differential conductance of 
both dots is the same, but the situation changes when a 
becomes smaller than 1, a < 1. Interestingly, the con- 
ductance of the dot weakly coupled to the leads (in our 
case the dot QD1) is larger than the conductance of the 
dot strongly coupled to the leads. For a sufficiently small 
value of a, the differential conductance of the dot QD2 
appears as a broad background, whereas the differential 
conductance of the dot weakly coupled to the leads is 
then very narrow. This behavior follows the features of 
LDOS of the dots QD1 and QD2 discussed above (Figg). 
Apart from this, we note that the total differential con- 
ductance diminishes as a decreases, and its line width 
also shrinks. This behavior indicates on the suppression 
of the effective Kondo temperature as a decreases. Such 
a behavior stems from the fact that the rate of tunneling 
events leading to the Kondo resonance decreases since 
the dot QD1 becomes detached from the leads as a de- 
creases. This is also in agreement with our predictions 
on the a dependence of the Kondo temperature, derived 
in Sections IIIII and llVl (see also Fig|2]). Finally, when one 
of the dot is totally disconnected from the leads (a = 0), 
the Kondo temperature vanishes, and no Kondo effect 
appears. 

When €\ ^ e-2 (eg. when the level slitting due to renor- 
malization is not compensated by external gate voltages) , 
the Kondo peak in differential conductance becomes split 
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FIG. 6: Differential conductance calculated for indicated 
value of a in the case when the level slitting due to renormal- 
ization is not compensated by external gate voltages. The po- 
sitions of the dots' levels have been estimated self-consistently 
using Eq. (J5j . The other parameters as in Fig[4] 



and the two components arc shifted from the Fermi level 
and have rather low intensity, as shown in Fig|51 This 
suppression of the Kondo anomaly resembles similar be- 
havior in the case of the spin Kondo effect. 

The presence of Kondo peaks also depends on the cou- 
pling strength of the dots to the leads. Above we assumed 
relatively strong coupling for both dots, with some asym- 
metry of this coupling described by the parameter a. We 
have also examined the case when both dots are weakly 
coupled to the leads. There is no Kondo effect in such 
a case, which is consistent with the recent experimental 
observations^. 



B. Numerical results for the case t 7^ 

Now, we consider the situation when direct hopping 
between the dots is allowed. Figures [7] and |8] show LDOS 
for the dots QD1 and QD2. When both dots are equally 
coupled to the leads (a = 1), a double peak structure 
emerges in the LDOS of the dots QD1 and QD2, and 
the LDOS is the same for both dots. The two peaks are 
centered at e = ±2i. This comes from the fact that when 
t 7^ 0, the dots' states hybridize into two molecular-like 
states with eigenenergies e± = eo ± t. These new states 
are then involved in the Kondo phenomenon. If initially 
an electron occupies the level e_, then it can tunnel into 
the Fermi level of a given lead and simultaneously an- 
other electron having energy +2t tunnels onto the level 
e_|_. Coherent superposition of many such events results 
in the Kondo peak at the energy e = 2t. In the same way 
one may explain the presence of the Kondo peak at the 
energy e = — 2t. 

However, the situation changes for a =/= 1. Apart from 
the two peaks located at ±2i, one finds an additional 
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FIG. 7: Local density of states for the dot QD1, calculated for 
indicated values of a and for t = — 0.1T. The other parameters 
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peak in the LDOS of the dot QD1, which is located at 
the Fermi level. However, instead of the peak, a dip in 
the LDOS of the dot QD2 appears at the Fermi level. 
Possible explanation of this behavior relies on the tran- 
sitions/tunneling events which do not induce electron 
exchange between the molecular states e+ and e_, but 
rather between original bare dot levels. The appearance 
of the dip (at the Fermi level) can be explained using 
the arguments from Section[V]A. One can notice that for 
q = 1 (symmetric couplings) there is no peak (or dip) 
in the LDOS at the Fermi level and LDOS for QD1 and 
QD2 are equal. For asymmetric couplings, i.e., a =/= 1, 
the amplitude of the LDOS at the Fermi level for the 
dot weakly coupled to the leads exceeds that for the dot 
strongly coupled to the leads. As a result, the amplitude 
of the LDOS at the Fermi level for the dot strongly cou- 
pled to the leads decreases (as a decreases) and the dip 
structure occurs. 

Differential conductance for several values of the asym- 
metry parameter a is shown in Fig[3] When the dots are 
connected in the T-shape geometry, a = 0, one finds two 
maxima centered at eV = ±2t, and one dip at eV = 0. 
Suppression of the conductance at eV = is a result of 
destructive quantum interference 31 . When the coupling 
to the dot QD1 is turned on, then the dip structure dis- 
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appears in the differential conductance. Instead of dip, 
one finds the third peak centered at eV = 0. For a fully 
symmetric system, a = 1, this peak vanishes and only 
the satellite maxima are present. It is also worth to note 
that the conductance increases with increasing a. 



ages, so one could reach the full Kondo anomaly. Such a 
compensation, however, is not possible when the direct 
tunneling between the dots is strong. 



VI. SUMMARY AND CONCLUSIONS 



We have considered the orbital Kondo effect in a spin- 
less system of two single-level quantum dots connected to 
electron reservoirs. Various techniques have been used to 
describe basic features of the Kondo physics. First, we 
used the scaling technique to evaluate the level renor- 
malization and the corresponding Kondo temperature. 
Then, we used the slave boson technique to calculate lo- 
cal density of states and linear conductance. To find non- 
linear conductance we used the nonequillibrium Green 
function method. 

The numerical results show that transport character- 
istics reveal typical Kondo phenomenon, similar to that 
observed in a single quantum dot with spin degenerate 
discrete level, coupled to external ferromagnetic leads. 
In the case considered, the splitting due to level renor- 
malization could be compensated by external gate volt- 
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Appendix: Green's functions 

Here we show explicit form of the derived dots' Green 
functions G- L j for i,j = 1,2; 

I 



Gu ~ M 



1 + — {njAj - n a t) 



(A.l) 



M {W 



t 7 / , ~\ 

1 + — (n^ - ra ?i tj 



(A.2) 



where 



M = ^11^22 - ^12^21 



Qu = e - e, - £f + — ( ai A- + j-t) 
U , 

Mil = t ~ {li A i + Oqt) , 



7l = S 22 + 



72 = ^22 + S n . 



The self-energies are defined in the following way 



= £ l^| 2 |(^ - ^ + Ae)F(")(e ka ) 



ka 



Ai = e - ei - U - - - Y, e 22 - S 



11 ^11 -^22 ^22i 



=^|^| 2 ^(e-e ka -Ae)F(")(e kQ ) 



ka 



A 2 = e - e 2 - U - £ 



(o) v c 

22 ^22 ^11 — ^11) 



t — t + + ^11, 



A 



a l — ^22 + ^11 + ^22 1 



a 2 — ^11 + ^ii + ^22, 



ka 



iv^l 2 



e + e ka - ei - e 2 - U 



\ - i T ,ai2 *( e - £ ka)(e - e ka + Ae) - 2t 2 M 



k Q 



/(«) 



a|2 *(e - £k Q )(e - e kQ - Ae) - 2t 2 , n 



ikl 



ka 



A 



£ko 



r 



Efl=E jiffli /n(etJ 



ka 



e - £ko 



with A = (e- e kQ ) [(e - e k « + Ae) (e- e ka - Ae) -At 2 } , Ae = 
ei -e 2 , Fi n \e ka ) = / a (e ka ) for n = J and Fi" } (e kQ ) = 1 
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for n — 0. The self energy E^ - is the self energy of the 
noninteracting system, i.e., U = in the Hamiltonian 

(I)- 

Assuming that e\ = &2 = £o the self-energies including 
Fermi distribution function can be calculated analytically 
and expressed by means of digamma function. 

Although, presented here EOM approach gives qualita- 
tively good (physical) results around Tk , it breaks down 
(both qualitatively and quantitatively) at lower temper- 
atures (especially at T <C Tk)- We must point out that 
the drawbacks in EOM method are due to the logarith- 
mic divergence of the di gam ma function. Detailed anal- 
ysis can be found in Ref. [33l - l36| . Here, we only list them 
briefly. Specifically, the divergence of the digamma func- 
tion leads to wrong behavior of the linear conductance 
(density of states at the Fermi level Slave-boson cal- 
culations (which are exact at T=0K) show that the con- 
ductance saturates as the dot's level is decreased (see 



Fig©. This is no longer true for EOM approaches, where 
below certain dot's level position the conductance starts 
to decrease 3 ^. The EOM methods, with different de- 
coupling scheme a 4 ' 33 ' 36 , also do not conserve complete- 
ness relation as well as not satisfy Friedel sum rule 3 ^. 
However, the decoupling schemes considered in our work 
somehow give surprisingly good dependence of the dot's 
occupation numbers as a function of the impurity level 
position (similar to those obtained by numerical renor- 
malization group method^). Taking all above into ac- 
count one should remember that presented here EOM 
method leads to incorrect predictions at low tempera- 
tures and is limited to T > Tk as we stated in Section 
fVl However, the EOM methods have an advantage on 
the other techniques used to investigate Kondo problem. 
Specifically, it enables to explore the non-equilibrium 
phenomena present in QDs systems biased by a finite 
voltage difference attached to the external leads. 
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